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Abstract 

This paper extends the calculation of quantum corrections to the cosmo- 
logical correlation (CO) which has been done by Weinberg for a loop of 
minimally coupled scalars, to other types of matter loops and a general and 
realistic potential. It is shown here that departures from scale invariance are 
never large even when Dirac, vector, and conformal scalar fields are present 
during inflation. No fine tuning is needed, in the sense that effective masses 
or coupling constants can have arbitrary values. Even when the mass is as 
large as Mpi, the one- loop result is still naturally smaller than the classical 
one. Thus, the results show that the existence of these fields during inflation 
may not be ruled out and is consistent with natural reheating. 
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I. INTRODUCTION 



If we would like to understand how electrons and photons arise during re- 
heating, it is natural to look back and ask what happened at the time during 
inflation. The gaussian and nearly scale invariant spectrum predicted by the 
scalar field dominated universe theory agrees very well with current obser- 
vations. It is therefore widely believed that the quantum fluctuations of a 
scalar field during inflation seeded the large scale structure of the universe 
we observe today. However, in order to understand where other matter such 
as fermions and photons observed today come from, the inflatorH would 
have to couple with these flelds during inflation. There is no reason why 
there must be only scalar flelds and gravity but nothing else during infla- 
tion. Other matter would not have arisen during reheating if there were 
only a scalar field that coupled exclusively to itself and gravity. If other 
matter such as Dirac, vector, and conformal scalar fields during inflation do 
not give anything far larger than the observed valued and do not break the 
scale invariant of the spectrum in the correlation functions, the existence of 
these flelds during inflation may not be ruled out. 

So far, the quantum theory of cosmological fiuctuations is considered up 
to the quadratic term in the action [1]. Recently, non-gaussian terms in 
the scalar field(s) have been calculated classically [7,8]. The quantum effect 
to arbitrary order in cosmological fluctuations has been more recently for- 
mulated by Weinberg [2]. With a sample massless minimally coupled scalar 
loop calculation, Weinberg's result shows that the momentum dependence 
goes as q~^lnq, with an additional suppression of GH^ when compared to 
the classical result. Is this true for other kinds of matter such as Dirac, 
vector, and conformal scalar fields as well? In fact, the unbroken symmetry 
matter becomes non-negligible when we go beyond the quadratic term in the 
action in cosmological perturbation theories. It is therefore of great interest 
to investigate how the higher-order corrections to the bilinear correlation 
function {COloops depend on momentum q when gravitational fiuctuations 

^Although we do not actually know what the inflaton is and is not, here we treat 
any field that has unbroken symmetries, i.e (x) = (V') = (^(j) = as other matter. A 
scalar field ip that has a non-zero expectation value is considered as an infiaton, as in 
conventional belief. Therefore in the quantum theory of cosmological fiuctuation during 
inflation considered here, the inflaton ip, gravity g^^, and other matter fields ^^'^ 
expanded as </p = <^ + S<fi,g^^ = g^^ + Sg^^,x = + (5x, V" = + 5^,Af, = + (5y4^, 
respectively. 

^The classical correlation function is (CC) = , 
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interact with general matter other than scalar fields. Will the result go ap- 
proximately as as in the scalar case? 

If other fields, other than inflaton, acquire a VEV, it is not necessary 
that the result will not affect the scale invariance of the curvature perturba- 
tion. For example, it is considered in the literature that Dirac fields cannot 
give a scale invariant primordial spectrum of density perturbations because 
their momentum dependence {(^Q is far from Q'~^[6] and vector fields can 
generate a scale invariant spectrum only in a special kind of mass [13]. How- 
ever, the Dirac and vector fields in [6] and [13] respectively were considered 
only classically. In general, Dirac and vector fields in an expanding universe 
only exist as quantum fields with zero expectation value {ip) = {A^) = and 
what we observe is the density correlation function related to (CC)j not the 
product of the fields {StpStp), {5Ai5Aj). Therefore, we have to learn how to 
quantize such fields at higher order in cosmological fiuctuations. In this pa- 
per I calculate the quantum effect to the observable {(^Q due to other matter 
loops. We use the in-in formalism [2, 4], appropriate for calculating expec- 
tation values, rather than the S-Matrix in time dependent backgrounds. I 
mainly follow the calculation of Weinberg for a loop of minimally coupled 
scalars in [2] and extend it to the loops of (massive) Dirac, vector, and con- 
formal scalar fields. We investigate how the ( correlation function depends 
on its momentum and whether it can be large. 

It is also important to investigate the momentum dependence when 
quantum corrections are applied to two-point correlation functions. If the 
momentum dependence of the loop spectrum goes as such that n is 
far greater than 3, this will produce a larger spectrum than the classical 
value outside the horizon when q ^ 0. The existence of Dirac and vector 
fields during inflation can be easily ruled out if the spectrum is far from the 
scale invariant; therefore, those fields cannot be the candidate for the origin 
of structure. We have shown that this is not the case. We always obtain 
nearly scale invariant spectrums even when Dirac, vector, and conformal 
scalar fields are included. 

In section II we explain why quantum effects could be large. In sec- 
tion III we present the aspects of the in-in formalism and the un-equal time 
(anti) commutators of all fluctuations that are needed for our present pur- 
poses. Section IV and V introduce a class of theories with a single inflaton 
field, plus an additional Dirac field with gravitational interactions. A (mas- 
sive) fermion loop correlation function of {(^Q is also calculated in these 
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sections. In sections VI, V II, and VIII we follow the same approach as in 
section IV and V except we replace the Dirac field with a vector field or con- 
formal scalar field respectively. In section IX we summarize all the results of 
this paper and explain why the departure from scale invariance is still small 
even in a more general potential V{(p) V((p,'ip'ip, A^A^^). All results show 
that for all theories and matter with a general potential V{Lp,'ip'ip, A^A^^), 
the quantum correlation functions are never much larger than the classi- 
cal (observed) value and are nearly scale invariant. In appendix we derive 
gravitational and general matter fluctuations used in the loop calculations 
of section IV- VIII to cubic order for the general reader. 



II. PROBLEMS 

There are some simple arguments that lead us to believe that quantum ef- 
fects might contribute to the spectrum in the order of the observed values 
without getting suppressed by an additional factor of G. This could happen 
when matter couples with the inflaton and gives a vertex in the order of ip. 
The fact that -jj^ is not small raises the question whether loop effect could 
be large. 



To clarify the problem, we take an example of the interaction of a fermion 
ip, inflaton cp, and graviton g^^i, via 

Cint = y/^'P'ipIp (1) 

In cosmological fluctuation, we normally expand the fields as 

g^_iu = g^iu + 5g^y,<p = 1^ + 5Lp,tp = ^ + Sil) (2) 

In general, fermions interact with fluctuations of both gravity 5g^u and scalar 
field 5ip and thus affect the conserved quantity Q. However, we can choose a 
gauge in which the inflaton does not fluctuate {6ip = gauge [7]) so that (" 
is purely gravity in this gauge. Therefore, one of the interactions in eq. ([T]), 
for instance the trilinear interaction, has the form 

H(^^^{t) = y"dW(t)^(t)VS(x,t)V(x,t)F[C(x,t)] (3) 

where F[Q is some function of C, depending on the details after the expansion 
of the metric. Let us choose F[C,] = C, for simplicity. Now we calculate the 
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quantum contribution of (CC)0 

c(x,t)c(x',t)) = - r r dh(\Hi, [i?2,c(x,t)c(x',t)ii) (4) 

' J —oo J —oo ^ '- '- J J / 

where Hi = Hj{ti), H2 = Hj{t2), and Hj is the interaction part of the 
Hamiltonian (the part that is of third or higher order in fluctuations) in the 
interaction picture. By solving Dirac field equation in an inflating univers^, 
the fermion pair ^pip goes as at late time. Therefore, the factor 
cancels with ^/—g for each Hamiltonian. With zero factors of a{t), the result 
of two time integrals grows as (Ina)^. But (Ina)^ is not as significant as ip 
in producing an appreciable effect in the interaction eq. ([3]). Since there are 
a total of four (s on the RHS of eq. ([!]), we get the factor |Cq |^ — ^^^^2^6 ^ ■ 
The two time integrals become 

where r is the conformal time r = Jj°° ^F)' Since (f does not change very 
much during inflation, it can be approximated as (p{ti) ~ (^(^2) — '^(^g) 
at the time of horizon exit. Thus, does not enter into the time integral. 
By collecting the factors of (p and SvrG, we can then approximate the 
correlation function as 

iCQioop j^2^^ (6) 

where e = — is a slow roll parameter and Cq is the momentum dependence 
that depends on the results of time and momentum integrals and the details 
of the propagator to the momentum p and p' . It is important to note that 
unlike the term H, the term (p, which could arise via Yukawa coupling, does 
not give a spectrum that is suppressed by a factor of GH^. Therefore, if an 
unperturbed inflaton amplitude is as large as Mpi, we have the correlation 
function 

which is in the order of the classical result. 



*See the full formula in the next section or in ref[2] 
^This is shown explicitly in fermion sections 
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The large vertices in the order of the Mpi raise the possibiUty that quan- 
tum effects arc not suppressed by a factor of G as was previously believed. 
Therefore the large vertices could contribute to the loop spectrum in the or- 
der of the classical value. However, such realistic <f{tq) ~ Mpi coupling can 
only happen in massive, but not massless, matter fields at one-loop level. 
The reason is that the inflaton fluctuates around a non-zero background that 
always contributes to the non-derivative matter terms in the second order 
after field expansions, i.e., m^'i^ = <f^ip or |<^P^f = m^Af. The massive 
case is more general because it allows the possibilities of interactions with 
other broken symmetry fields such as the inflaton and hence gives large ver- 
tices in the order of Mpi. Although the argument above is valid, we still 
need to find out what Cq is through actual calculations because Cq is also a 
function of mass that arises through the massive propagators. We have to 
investigate whether this Cq{m) gives other kind of suppression. 



Since the Hamiltonian that governs fluctuations of fields is explicitly time 
dependent, we need an in-in formalism[4]. For the purpose of calculating the 
in-in correlation function to arbitrary order, it is convenient to use Wein- 
berg's formula[2] 



The expectation value on the RHS of the equation above is annihilated in 
the free field vacuum by annihilation operators and on the LHS it is in inter- 
acting vacuum. We will use the equation above to calculate loop correlation 
functions. 

We see from Weinberg's formula above that quantization requires un- 
equal time (anti) commutators. The quantity C(x, t), the scalar field x(x, t), 
the Dirac field 'i/'(x,t), and the vector field Ai{-x.,t) can be Fourier decom- 



III. CORRELATION FUNCTION FORMULA 





x([iJ,(ii), [Hi{t2),...[Hi{tN),Q'm-]])^ (8) 
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posed as 
C(x,t) 
X(x,t) 

A(x,t) 



J d\[e'^-McDC,it) + e-*'i-a*(q)C*W] (9) 
J d'q[e'^-M^)xS) + e-'^-"a*(q)x;(t)] (10) 
f d'qY, p'»-"«(q, s)t^q,.(i) + e-'^^-^PH^, s)V^At)] (11) 

s 

J d'qY^ p'i-e,(g,A)a(q,A)^,(0 + e-*'i-e|(g,A)a*(q,AM*((l|) 



where s = it^ stands for the spin of fermion, A = 1, 2 is a hclicity index of 
a photon, A = 1, 2, 3 is a hehcity index of a massive vector field. ei{q, A) is 
a polarization vector. The creation and annihilation operators satisfy usual 
(anti) commutation relations as 



«(q),«*(q')J = [«(q),a*(q')J ='^'(q-q) 

a(q,A),a*(q',A')] = 5AA"J'(q " q') 
{a(q,s),at(q',s')} = {/3(q,5),/3t(q',s')}=5,,,^3(q-q') 
{a(q,s),a(q',s')} = {/3(q, s), /3(q', s')} = 
a(q),a(q') = a(q),a(q') = Q!(q,A),a(q',A) 



(13) 

(14) 
(15) 
(16) 
(17) 



and Cq{t), X(/(0' Uq^sit)^ ^,s(^) ^-nd ylq(t) satisfy interaction picture free field 
equations in an inflating universe. With the (anti) commutation relations 
above, the unequal time (anti) commutators of all matter are 

C(xi,ti),C(x2,t2)] = 2ij dVP'("^-"^^/m(Cp(ii)C;(i2)) (18) 
X(xi,ii),x(x2,i2)] =2iJ d'pe'^<^'-^'hm(^Xp{h)Xp{i^)) (19) 

{V'(xi,ii),VS(x2,t2)} = J dVP'^"^~"'^l^(t^p,.(il)^^p,5(t2)+Vlp,.(tl)F-p,s(t2 

(20) 
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and 



^i(xi,ti),A,(x2,t2)l = 2i / d^pe'P<^'-^^^U^jIm(Ap{h)Al{t2)) (21 



where Ilij is the polarization factor in which 



A=l 

is time independent for m = 0. 



PiPj 



(22) 



In this section, we calculate a loop power spectrum of general matter 
valid for scalar, fermion and gauge fields. To see this in more detail, let us 
consider a general complex bosonic or fermionic field ^' with the interactioii^l 



Hi{t)= j d^xV{t)^*{x,t)^{x,t)C{x,t) 



(23) 



where we omit the spinor and vector indices in this section, and V{t) is any 
time dependent vertex function as a consequence of an expanding universe. 
The loop correction to the correlation function when = 2 is 



Q{t) 



dU 



t2 



dU 



H2,Q 



(24) 



If Q is the product of the gravitational field ^ treated as the external legs, the 
vacuum expectation values of matter, either bosons or fermions, that circu- 
late inside the loop on the RHS of eq. (p^ can be evaluated independently 
from Q because they have different types of creation and annihilation oper- 
ators. We can write a real C, field and a complex ^ field in the interaction 
picture in terms of creation and annihilation operators as 



(25) 



C(x,t)= / d^p e^P-apCp(t) + e-^P-XC;(t) 



M/(x,t) = I d^pY, (e'P-''«p,AXp,A(t) + e-^'^-^pl^W^^xit)) (26) 



We can easily generalize to realistic interactions (such as the terms containing field 
derivative) later in the next section, once we obtain a general one-loop formula valid for 
any matter in this section. 
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and 



^*(x,t) = / d'pY, (e~^P-"<A^;,A(i) + e'''"'Pp,xWlxit)) (27) 

where Op satisfies the commutation relations since is a boson but ctp^x and 
Pp^x satisfy the (anti) commutation relations for (fermionic) bosonic matter 
loops. A is either the spin index for the fermion or the helicity for the gauge 
field. Xp^x{t) and Wp^xit) satisfy the matter free field equation in an inflat- 
ing universe. 

Let us now evaluate the RHS of eq. ([2 



H2,Q 



d'x2V2 



^2^2C2,Q = / d'a;2F2^'2^2(C2Q - QC2) (28) 



Hi, 



H2,Q 



(fxid^X2ViV2^ \^\^lCl,'^l^2C2Q 

vl/*M/iCi,Qvi/*M/2C2l) 

J /o 

j j <fxid'x2ViV2(J^n'^in'^2{ClC2Q - C1QC2) 
+ ^'*2^2*l^l(QC2Cl - C2QC1))) (29) 



Since is independent of C and Q, the vacuum expectation value can be 
evaluated independently. 

For the C, part, 
Hence, 

(CiC2Q)o = (QC2Ci)S (31) 

Similarly, 

(CiQC2)o = 2 j d3A;d3fe'e*'^-(— )+*'^'-(-^"-')a(ti)C-fc'(t)a(t)Cfc'(i2) (32) 



Hence, 



(CiQC2)o = (C2QC: 



(33) 
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For the matter field, only ^'i could pair with ^'2 at different times. To see 
this in more detail, we write the field operator ^ in terms of creation and 
annihilation operators similar to that in eq. (j26p . Hence, these products of 
the fields can be written in momentum space as 

X 5]Xp,A(^l)x;,,(^2)H^p^A'(^2)T^*^A'(^l) (34) 

A, A' 

Since 

= (^*^r2^*^r^)* (35) 

Hence, the correlation function is 

J d3xe*'i-(^-^')(|Q(t)^ = - J' dt2 p dti j Sx j Sx2 j Sxi 

x2Re({^\^i-^l^2)^{{CiC2Q)» - (CiQC2)o)) (36) 



Note that many integrals over the momenta k and k' can be eliminated via 
the space integrations that produce delta functions, i.e.. 



d^x (2^)3(53(q-k) (37) 

j d^xi ^ (27r)353(q + p + p') (38) 
j d^X2 ^ (2^)35=^(k'-p-p') (39) 



We therefore have the formula 



J (i3xe*'i-('^-"')(c(x,t)C(x',t))^^^^ = -4(27r)9y d'pd^S^iq+p + p') 

X / dt2V2 r dtiViRe(zM) (40) 

J —oo J—oo ^ ' 

where 

2 = Q{tiK:{t){Q{t2)c;{t) - c,(t)C(t2)) (41) 

and 

M = J2Xp,x{ti)Xlx{t2)W,,,y{t2)W;,^y{h) (42) 



A,A' 
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Eq. (j42p is a formula for a general matter field loop. 



For a real scalar field x(x, t) = t), we have = Xp and Wp = Xp- 
Hence, eq. ([12]) becomes 



-Mx = 2Xp(ii)x;(i2)x;'(i2)Xp'(ii) (43) 
Note that we have an additional factor of 2 for any real field. 

For the charged scalar field t) / x*(x, t), we still have Xp = Xp and 
VFp = Xp- Hence, eq. (|i2]) becomes 

-Mx = Xp(ii)x;(i2)x;'(i2)Xp'(ii) (44) 
which only differs from the real scalar field by a factor of 2. 

For a fermionic field which is always complex, we have Xp^\ = Up^s and 
Wp^x = Vp^s- eq. (02]) becomes 

Yl UpAh)u;At2)Vp,,s{t2)v;,Ah) (45) 

r,s=l,2 

For a real vector field j4j(x, t) = j4*(x, f), we have Xp^x = Ap^x and Wp^x = 
^*;^. Hence, eq. (0^ becomes 

Ma = 2Y,AAtM*p,xit2)A;,^x'it2)Ap>,x'iti) (46) 

A, A' 

Note that the fermion and gauge fields require the summation over spin 
and helicity at different times, especially in the massive theories. In the 
subsequent sections we will use the formula for the general matter field loop 
shown above to calculate the power spectrums for more realistic interactions 
between various kinds of matter and gravity. 

IV. FERMION LOOP, INFLATON, AND GRAVITY 

In the known theory of cosmological fluctuation, only an inflaton and gravity 
are considered. Here we consider additional fermion as 



C — Cg + Cip + Cf 



1 , 

9 



2 



SvrG 



+^j''{Vo,^p) - (PaV^)7"V' + 2mVSV (47) 
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where /i, v, ... are the space time indices, and a, (5, ... are the Lorentz indices 
raised and lowered by the vierbein V^. To deal with fermion, we need 
the tetrad formalism[10]. The metric in any general non-inertial coordinate 
system is related to the vierbein by 

9A^) = V^{x)VS{x)r^ap (48) 
The covariant derivative to fermionic field due to the gravity interaction is 

-Da = V^d, + \aP^Vl^V>iV,,.,, (49) 

where o"^'*' = |[7^,7'^]- To find out what the time dependent propagators 
are, it is necessary to solve interaction free field equations in an infiating 
universe. For a fermion, it is 

a-i-^('at7V') + — V' + "^V' = (50) 
dt\ ) a 

We can re-scale the field ^ = a'^S and work with the conformal time. We 
therefore have a simpler Dirac equation as 

7° 5' + -f'diS + maS = (51) 

where S' denotes the conformal time derivative. Since the background is 
spatially translation invariant, the solution can be written in mode function 

as 



V'(x,t) = a~i(t)S'(x,t) 



2 (t) J „ L 



32) 



a2{t) 

where Uq,s(t) and fq,s(i) satisfy 

7°^q,s + n\iU(i,s + imauq^s = (53) 

and 



7°^q,s - iY<liV(i,s + imavq^s = (54) 



For Ccj(^)) it satisfies Mukhanov's equation [1] as0 



C, + |(lna3(t)e)4+^C, = (55) 



^We work with interaction picture free field equations here to obtain time dependent 
propagators via their solutions. When loop effect is included, Mukhanov's equation is 
modified by varying the loop quantum effective action with respect to 
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We can expand matter and gravity fluctuations in the actions of (j47p to 
arbitrary order in cosmological fluctuation. It is also convenient to write 
down the gravitational and all matters actions in ADM form [5] and solve for 
and Ni in the constraint equations. The cubic term and higher order terms 
are time dependent vertices that are needed to calculate loop diagrams. The 
direct expansion of matter and gravitational fluctuations to higher order 
are complicated. However, as shown in the appendix, many terms are not 
necessary since they are cancelled via field equations and are removed by 
the field redefinition of C. Therefore, the important terms of the trilinear 
interaction of any general matter are 

^CMM(i) = - J d^xeHa^iT^'^ + a^T')V-^C (56) 
The fermion energy momentum tensor in the presence of gravity is 

Tf" = --{'lijf'Yi' + Vi^i^ - i^i^i^'^ - i^Y'^^) (57) 

Therefore, we have the time component of the energy momentum tensor to 
quadratic order as 

= {-i)[jj^ii)-'4)^^ (58) 

where -0 = ■0^/?, /3 = Z7°, (7*^)^ = —1, and = — 5o = — ^. Similarly, the 
spatial component of the energy momentum tensor is 

a^Tf = i(v^y(9,V')-(a,v^)7V) (59) 

where 5* = ^ and 7^ = Va^°'. Therefore, 7*^ = Vq'7*^ = 7*^ and 7* = 

1^7"^ = "J . These give us the Hamiltonian interaction of fermion and 
gravity to the cubic order as 



H^^^it) = -J d'xeHa' 



^(Vif (9^^)-(5,v;)7V)]v-2c (60) 



+ 

a 

Note that the interaction above is real and also valid for massive fermions. 
The cubic order in the interaction above can be further simplified with 

ijj^ip-i)j°ijj + —Yi;-^j'—-2m'ipip = (61) 
a a 
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Hence, the trilinear interaction Hamiltonian is 



= -2 f d^xeHa^l^'ip-f^i; -ijjj^i; + m^Pi?jV-^C (62) 



We see that Dirac's equations in an expanding universe (|53|) and (|54|) at 
m = are the same as those in Minkowski space except that the physical 
time t is replaced with the conformal time r. So we can expect the plane 
wave solutions for Uq^s(t) and Vq^s(t) to be 

^q,.(i) = n°^^,e-"^-,v^,s{t) = ^.e*"^^ (63) 

where u^ g and v^ g stand for constant coefficients outside the horizon. These 
coefficients can be determined by matching the solutions at deep inside the 
horizon with the flat space solutions. Deep inside the horizon, the field 
does not feel the effect of the expansion of the universe. Therefore, the 
normalization factor can be chosen in the same way as in Minkowski spac^ 



where q = q'^ = \/q^- We see that the momentum dependence q of the 
expectation value of the fermion and anti-fermion pair (tpip) is far from the 
scale invariant spectrum q~^. However, this does not rule out that fermions 
could not seed the large scale structure of the universe observed today. The 
reason is that we never observe the product of either scalar or fermionic fields 
but rather the product of temperature (^^) or density (^^) fluctuations, 
which are related to the conserved quantities (CO- Since fermions interact 
with the gravitational fluctuation, we therefore calculate how fermions affect 
at the loop quantum level. We now calculate the one-loop graph with two 
vertices of the two point function. Owing to the interaction of gravity and 
fermionic fields, the quantum corrections to the C spectrum are 

C(x,t)C(x',t))^^^^ = -J' dt2p dh(^[H,,[H2,a^,t)C{^',t)]')(65) 

For simplicity, we first start with massless fermions. We can use the formula 
(|4U|) and (j45|) derived in the previous section. To match with the interaction 



We emphasize that this formula is only valid for massless fermion. The situation is 
entirely different for massive fermion, as wiU be shown in the next section. 
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Hamiltonian in eq. ()62p . we replace the interaction in eq. (|23p with 

^ (iPj^^^P-(^ipWip) (66) 
\ a a ) 

Uti,2) ^ -C,(ti,2)/g' (67) 

y(t) = 2eHa^{t) (68) 
Hence, Z in eq. (HTI) changes to 

^ - ^(4(ii)C(t)(4(*2)C(t) - C,(t)C*(t2))) (69) 
and M.^ in eq. (|15]) changes to 

r,s 

The equation above shows the need to sum over spins at different times. 
Fortunately, massless fermions are conformally flat so we can still use the 
spin sum formula from flat space. As seen from Dirac's equation, the 
solutions of massless fermions are just plane waves with conformal time, 
'Up,s(7") = Up ,,e~^P'^, after re-scaling the field such that Up^s = CL~2{t)up^sit)- 
With the spin sum equation ([6^ . eq. ([70|) becomes 



= (p. - p[){p, - P;-)^(^J^^tr(f 7'^7^7^)e-(^+^')(---) (71) 

Since 

= 4(V°7?^'^ - V^'r/"^ + r/^^r/"^') (72) 

we have, 

= ,^ .g 4 4 (p-pO'(1 +p-p')e-^(P-^P'^(^^-^^^ (73) 
(27r)°a|a| 

Substituting into eq. (gO]) with y(t) = 2eHa^{t), we have 

y d=^xe^'i-(^-''')^C(x,t)C(x',t))^ =-16{2TTf J dVV<5^(q + P + p') 

(p_p')2(i+^.^') /■* dt2e2H2a2 r dheiHiaiRe(Ze-"^P+P'^^^'-^^A{7A) 

where ^ is the contribution of the ^ part in eq. (j69p . To calculate ^, we use 
the solution of the free field Mukhanov's equation in the interaction picture 

Q{t) = Ce-^'^^il + iqr) (75) 
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where 



Hence, 



2 ^ ^^GH\t,) 



|('o|4 



HiH2ala2 

During slow roll inflation, we approximate ei ~ 62 ~ ^{tq)- Integrating over 
conformal time ti, we get 

I d3xe*'i-(--')^C(x,t)C(x',t))^^^^ = -16(2vr)3|C°re2(t,) J d^P 

x6HcL + P + p'){p-p)Hl+p-p)Re r dT2 , \ (e-^-^-2_i) (73) 

J-00 g + p + 

where an upper limit t — > cx) or r — > means a time still during inflation but 
sufficiently late so that a{t) is many e-foldings larger than its value when ^ 
falls below H. Integrating over conformal time T2 gives 

Re f dT2 \ , {e-^'^^- - 1) = , ^ ^ „ (79) 
J -00 q+P+P 2q{q+p + p) 

Substituting eqs. (f76]l and (f79l) into eq. (f78l) . we have 



X / dVV'^'(q + p + p) ^^"^'^ (80) 

J q+p + p' 

Power counting shows that the results of the momentum integrals p will 
go as q^. If we use dimensional regularization to remove UV divergences, 
the finite part of (CC) for the massless fermion loop will go as q~^ In q. To 
determine the coefficient of the finite part of the momentum integral above, 
we will follow the calculation as done in [2] for the scalar case. Note that, 
in general, 

j d'p I dyS^ip + p' + cOf{p,p',q) = J d'pf{p,P = \p'\ = |p + qU)(81) 



and 



/POO /"l /'27r 

d^p= / p'^dp / d{cos9) / dip (82) 
Jq J-1 Jo 
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The conservation of momentum = |p + q| gives 



p'^ = + + 2pq cos 9pg 



(83) 



where 9pg is an angle between the vectors p and q. Since g is a fixed external 
momentum, we can choose in the z— direction. Hence, 9pq = 9 and 

p'dp' = pqd{cos 9) (84) 

With eqs. ([HI]), ([82]), and dH]), we have 

d-^p / d^p'6^{p + p' + q)f{p,p' ,q) = — / pdp p'dp f{p,p' ,q) {85) 

J Q Jo J\p-q\ 

Since 2p ■ p' = q^ — p^ — p'"^, eq. (j80p can be written as 

^3^g*q.(x-x')/^(x,t)C(x',i)^ 



_ 2{87rGH^{tg)f r27r ^ 
loop (27r)3g7 L g ^ 



where 



pdp 



/ , / {P-P') 



/\2 



\p+q\ 

p'dp 

\p-q\ q + p + p' 



1 + 



2 2/9 

2 — p — p^ 

2^?y 



(86) 



(87) 



With dimensional regularization, the UV divergence of the integral above 
for 5 = {) gives the pole term as 

27r 



-^J{q) /+^F(5) = q^e^^'"'F{6) 



where 



F{6) ^ ^ + Fi 





Therefore, in the limit 5 — > 0, 



-Jiq)=q'(^+Folnq + F^ 
q \ 



Frilnq + L 



(89) 



(90) 



where L is a divergent constant. To calculate the coefficient Fq, it is neces- 
sary to evaluate the integral ([87|) . The integral over p' give^ 

JiQ) = 2 y dp[p{15p'^ + 17m + Vy - ^(ll^' + 6g) 



'0 



4p(p + q){2p + <7)2 logb + q + p']) (91) 



'We use Mathematica for the integrals and six derivatives. 
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To eliminate the divergence in the momentum integral above, we differenti- 
ate J{q) in eq. (|9ip six times and then do the integration over p. This gives 
a finite result as 

where we take the limit q ^ after integrating over p. Hence, 

J{q) = qH-^^^Q + L) (93) 



5! 
or 

Fo = -- (94) 

Substituting J{q) back into eq. ()86p . we have the finite part of correlation 
function as 

/ /.e-.--.(c(x, .,((x'. ,),^^ . -^^m^ [.n , . C] ,.5, 

with C, an unknown constant. Notice that we have the same sign as that in 
the massless scalar loop because we do not have the time ordered product 
of fermion pairs in eq. ([5]). The opposite sign of fermion loops only arises 
in the in-out theory when we time order the product of fermion pairs in 
order to close the loop. Moreover, the result in eq. (195p is smaller than the 
classical result by a factor of GH^. 



V. MASSIVE FERMION 

The calculation is much more difficult for massive fermions. This is because 
the mode solution of a massive fermion at arbitrary wavelength during infia- 
tion is not a simple plane wave like in the massless case or in flat space. We 
cannot rely on the trace technology normally used for spinors in flat space 
since the spin sum at different times cannot be written in the compact form 
of 7 matrices. We therefore need to solve the massive Dirac equation in an 
expanding universe during inflation and perform the spin sum at different 
times by multiplying the matrices. 

We define Us,p(t) = [n+,p(r)S', n_,p(T)S']^ and U5,p(r) = [v+^p{t)S,V-^p{t 
where S are the two component eigenvectors of the helicity operators. We 
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use the Dirac representation of the gamma matrices. 

7° = (%) 

- (-i) ( ) m 

Therefore, eq. (j53|) gives 

p + i((7 • p^Uzp^p lb imau±^p = (98) 

We see that the equations above are first order coupled differential equations. 
To decouple them, we differentiate those equations one more time. With 
some algebra, we get two uncoupled second-order differential equations, 

n'4 p + {f + {maf ± i(am)')n±,p = (99) 

where {5 ■ pf = {pi + + pi) (^^^ 1 ) ~ ^-''2x2 = P^- Therefore, the 
equations above are solvable as 

ul^p+(f + X^{r^±ir)Y±^p = Q (100) 

The solutions are 

UFRw{V,S,t) = 



(101) 



n+,p X 5* \ ^ / u^^p X S 



u-,p X S J \ {p- ^)ufi,p X S 

I c+,j,V^Hj}\-pT) X S \ 
[ {p-a)c.,pV^Hj^\-pT)xS J 

where = | — ir, /i = ^ + ir, and r = ^. 

We choose the initial conditions so that the positive frequency mode 
solutions match with the flat space-time solutions deep inside the horizon 
[150 

E + m\h( IxS \ ^,j^Eit')dt' (102) 



^°Note that the asymptotic behavior of Hankel's function for large x = —pr is Hi, {x) 
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where E"^ = m? + S = ^ q ^ ^P^^-^ ~ ( 1 ) ^v"^"^ 

down. Note that S'^ S = 1. We can check that, in the massless hmit, 
""/'li = ^^^=Tf expected. 

To find the normahzation coefficients c±(p), we match the solutions of 
eqs. (jlUip with pU2|) by using the asymptotic property of Hankel's functions. 
Hence, 



IJT^V -LIT'- , 

c±,p = (103) 
2(27r)2 

To calculate we use the formula (f40l) and (l45]l derived in the previous 
section. To take account of the more realistic interaction Hamiltonian that 
arises after the expansion in eq. (j62p . we replace the interaction in eq. (|23p 
with 

^ (^VS7°V' - ';/'7°V' + ^V'V') (104) 

C,(ii,2) ^ -4(*i,2)/9' (105) 
y(t) = -2eHa^{t) (106) 

2^ is still the same as that in eq. (|69p . However, for massive fermion in 
eq. (P5]) becomes 

2 

= ^v^^7fl (107) 



l"'2 

where the as are 

a = Up,yUij,^p - Ufi^p'U^^p + imup_^p>Ufj_^p (108) 

a = a{m —m) (109) 
and Uf^^fi are the massive mode solutions of Dirac equations 

u,{x) = ^^e'fHW^x) (110) 



2(27r) 



M^) = -^^e—H}t'{x) (111) 
2(27r)2 
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for X = —pT. The factor arises due to the fixing of the coefficients c±^p 
of mode solutions deep inside the horizon with the flat space-time solutions. 
We can check that in the massless limit m = 0, u,, = Ur, = ui = — rc"*^"^ 

and a{t) = a{t) = ^^{p' - p)e-'(p+p')- . Hence, 

4 

M\m=0 = -^(^2^l{'i^ + P ■ P') 



1"2 



^ (p-p')^(l+^5-p)e~^(P+P')("^~"2) (112) 



(27r)6afa| 

which agrees with the result in massless fermion section in eq. (j73p . 
Substituting eqs. (j76p and (j77p back into the power spectrum formula eq. 
(liOl) with V{t) = -2eHa^{t), we have (CC) due to massive fermion loop as 

/d».e-<'-vc(.,,)c(x'.*)) = 

J \ / loop q 

X 



j d^pd^p'6\q + p + p')I (113) 



where 



/oo rt2 
dt2al{e-'^^^ - e*«^2) / dtiale''^^^ (114) 
-oo ./ — oo 

where is the result contributed by the fermionic part in eq. (jl07p and 
the exponential terms e^*^"^ in eq. (jll4p come from the parts. 

So far, the result in eq. (jl07p is exact. There is no approximation in- 
volved in eq. (|107p . However, the exact result involves the integrand of 
Hankel's functions and their time derivatives and this makes the integration 
quite challenging. Nevertheless, we get some idea that the time integrals 
will converge and go at most as (loga)^. However, in order to determine the 
momentum dependence of the power spectrum, we need to integrate over 
unequal times ^1,^2 aiid momentums p,p' associated with fermion fluctua- 
tions. We will first integrate over times and then over momentums. The di- 
rect calculation is complicated because it involves integrating over products 
of Hankel's functions with complex order. Since p is the running momentum 
from to 00 whereas q is the fixed external momentum associated with a 
conserved quantity ^g, it is helpful to divide the integrals over momentum 
p in eq. (jllSp as an integral when Ag < p < 00 and < p < Ag. The first 
integral when p — > cxd can be approximated as if the fermion is massless due 
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to its high momentum. So, the result at high momentum after dimensional 
regularization will be close to that in the massless case in eq. (|95p . The 
second integral, when p < Aq, indicates that the mass effect may become 
important in the result. Therefore, additional calculations are needed to 
determine the momentum dependence. 

The momentum p corresponds to the fermion fluctuation ipp and the mo- 
mentum q corresponds to the conserved quantity (g fluctuation, this implies 
that a massive fermion ipp exits the horizon before Cq exits the horizon if A is 
in the order of 1. In other words, the fermion fluctuation a^t/^piipp is frozen 
by the time the fluctuation (g crosses the horizon. Outside the horizon of 
Cg, p,p' are sufficiently small and Hankel's function can be approximated as 

111] 

which is valid for x <C 1 and /3 > 0. Hence, the mode solutions in eq. (jllOl) 
become 



Trr 



With this approximation, their (conformal) time derivatives are 

u° = -ul,uP = --ul (118) 



Hence, 



a° = -imul^,u°^^p (119) 



^*on ^*o-o m^|r(^)|'* /'101^ 

„*o-o (i^*o„o\* |1^(m)I i2ir~2ir /'100^ 
= ('^2 C^l) = ^ ^ ^ ' 

Notice that (T{p,p') and a'{p,p') are both time independent. The exponential 
factors e^'^'' arise when we fix the coefficients to match the solutions with 
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those inside the horizon. Hence, from eq. (|107p 



2m2|r(^)|4 
(27r)8afa-2 

3„— 3 -r-o 



Ml = : :v {2-{p-p'){p"-p-'-+c.c.] 



= ar'«2'-^pV (123) 

where the factors from the fermionic part will cancel with the factor 
\/—g in eq. (jll4p . Therefore, 



/oo nt2 
-oo J — oo 

Although we see from eq. (jl24p that the time integral is of the order of 
(loga)^, we need to evaluate this integral if we want to see the momentum 
dependence q for the correlation function of C,. 
From eq. (jl24p . we have 

= ^TTT^e r ^(e-?-2-e*^^^)ii;i(-z(?r2) (125) 

\tq) J-oo 

Using Ei{—ix) = ci{x) — isi{x), we have 

^" = -iSf^ f — sm(gr2).z(gT2) (126) 
With Mathematica, 

/ sm(gT2)s^(^r2) = — - (127) 

J-oo 7-2 



Therefore, 



T° = ^ (128) 



Notice that the result of the time integral above is (/-independent. Hence, 

rp+q 

'\p-g\' 



/ p'dp'J^;,^, = G"{p + q)- G°(|p -q\)= 2q—- + 0{q^) (129) 

J\p-q\ CIP 
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Note that we only keep the leading order in q in the last equation. Hence, 

— / dpp dppi = — — (130) 

Q Jo J\p-q\ 3(27r)«F^ 

Substituting eqs. ()130p into ()113p . we have the C correlation function due 
to a massive fermion loop as 

Using |r(/i)|2 = |r(i ± ir)|2 = we have 

which goes as at low momentum. 



It should be noted that the power spectrum will not be zero when we 
take m = 0. Eq. (jl32p for massive fermions is the result when the fermion 
pairs 'ippiipp exit the horizon before or at the same time as C,q exits the hori- 
zon {p < q). We keep only the most dominant mode solutions for massive 
fermion after horizon exit. For the massless fermion case, the solution is 
simple enough so that we can do the integration exactly without any approx- 
imation. The integrand contributed by the massive fermion !F° = a\a^M.^ 
becomes frozen after horizon exit. The negative power of (— r) that arises 
in the time integrals j J dt2dti of the massive fermion, but not the mass- 
less fermion, is more important than the exponential function when r ^ 0. 
Therefore, massive fermion loops can contribute the (loga)^ factor because 
the interaction goes as a^. In comparison, logo does not arise in the mass- 
less fermion case because the interaction goes as rather than a°. 

As mentioned in section II, a fermion mass could arise from the non-zero 
vacuum expectation value of an inflaton field in a flat potential. Therefore, 
the effective fermion mass during inflation could be as large as Mpi. We 
know that in order to generate all matter observed today, the inflaton ip 
must couple to matter such as fermions sometime during inflation. Since we 
work in the gauge where an inflaton does not fluctuate 5(p = 0, the Yukawa 
coupling that can arise for the general inflaton potential does not change 
the result in eq. (|132p but only shifts the fermion mass to be 

m — > m + ifitq) (133) 
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However, during inflation, the fermion mass could be large because the non- 
zero expectation value of the scalar field ip{t) could be large. If the unper- 
turbed inflaton amplitude at the time of horizon exit is as large as Mpi, 
we have m ~ Mpi = ,^ Therefore, the power spectrum due to massive 

V SttG 

fermion loops is 

/ d'xe^^<---'Uc{^,t)C{^',t)) = ^^^""^l^l (134) 

J Wloop,m=Mpi 24g3cOsh2^ ^ ' 

We see that even when we include the large (p ~ Mpi coupling that seems 
to give the quantum effect that does not get suppressed by the factor of G, 
the result is suppressed by the factor cosh^ Mn2L instead. This happens be- 
cause the inflaton fluctuates around a non-zero background, implying that 
the massive fermion propagators are of the order of Mpi . The factor cosh ^ 
arises when we fix the mass dependent coefficients of the mode solution at 
late time. 

It should be mentioned here that the large mass term does not get 
suppressed to the quantity like a'^{'ipijj) — > tanh^. Terms like uu or vv 
approach constants with the mass dependent constant coefficient going as 
tanh which does not have large mass suppression. However, to close the 
fermion loop, the tree {ipip) is not the only quantity we need to calculate. 
We also must calculate (V'iV'iV'2V'2), which is the trace over the multiplied 
matrices such as s ^p' ,r{'ti)up,s(ti)up^s{t2)vp' ,r{t2) ■ Therefore, the bilinear 
viui gives a constant mass dependent coefficient of ^oGh ' resulting in a 
small one-loop result even when the fermion couples to the inflaton in the 
order of Mp;. 

Apart from ultraviolet divergences, no infrared divergence can arise due 
to the late time behavior. The reason for infrared safety comes from the 
fact that the function approaches a constant at low momentum. This is 
similar to viewing a and a in eq. (jl07p as scalars which approach constants 
after horizon exit. Provided that the integral over time is infrared safe, the 
integral over time only comes from the C correlator, whereas the fermionic 
part only contributes an factor that always cancels with the factor ^/—g 
in each interaction Hamiltonian. After all integrations, the power spectrum 
gives a momentum dependence. The result of the massive fermion case 
is valid at low momentum modes only where we need to cut off the momen- 
tum integral p to some value i.e., Aq, so that the approximation of small 
p,p' in the Hankel function (eq. (jllSp ) is still valid. This case means that. 
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if A is of order 1, the fermion momentum modes exit the horizon before 
the C momentum mode q crosses the horizon {p,p' < q = a{tq)H{tq)). For 
higher momentum modes p, the fermion behaves hke it is massless and is 
always suppressed by the factor G and a negative power of the Robertson 
Walker a as shown in massless fermion section. 



To investigate whether the quantum effect is truly small, more careful 
consideration is needed for the mass effect. The reason is that various mass 
dependent coefficients can arise when matching the general solution with 
that inside the horizon for the general graphs. However, a fermion has two 
components that are needed to form a pair with its conjugate. Bilincars like 

__ _ ni/\i9/ 7rm 

uu or vv, but not vu, contribute a constant factors like |r(/x)| (e h —e h ) 
which are tanh ^ — ^ 1 in the large mass limit. Bilinears like u^u or v^v 

cosh 

contribute a constant factor like ^^^-^ Mtt = 1, which is mass independent. 
Therefore, by considering this alone, a fermion has no exponential suppres- 
sion and seems to give a large quantum effect if a vertex is as large as Mpi . 
However, as shown in the detailed calculation here, this is not possible for 
the loop graph that has two external legs with two trilinear vertices because 
it requires a bilinear like vu instead. Bilinear terms like V'7V get suppressed 
at late time because 7* can only have a contraction with ^. Therefore, the 
result is suppressed by an additional negative power of a and its low mo- 
mentum outside the horizon. Interaction terms like 'ij^^^ip contribute both 
^?/;7*'0 and m'ip'ip factors via Dirac equation and its conjugate. Therefore, 
the maximum result of ijj^f^tjj cannot exceed the result of V'V'- The other 
powers of bilinear terms like (tljip)^ for n > 1 cannot couple to the mass 
dimension in the order of Mpi and can only give higher fermion loops (by 
dimension counting in the action). Hence we expect the ['ijjTp)^^^ interaction 
type to be suppressed by a negative power of a{t). 



VI. VECTOR LOOP, INFLATON, AND GRAVITY 

Ford considered a classical vector field driving inflation [12]. In this section, 

we consider a quantized vector field that affects the quantity and its corre- 
lation function through the interaction with gravitational fluctuations. The 
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action is 



SvrG 



(135) 



To determine the vector field propagator, we need to solve the interaction 
free field equation in an inflating universe. It is 







For 1/ = 0, this gives 



For I' = j, this gives 



An 



IQiAi 



Aj + HAj + ^Aj+m'^A 



_ qjQiAi 



■J 



+ iqj{Ao + HAo) 



(136) 



(137) 



(138) 



To eliminate the auxiliary field Aq, we apply to eq. (|136|) . We have 

d^ia^A'') = (139) 



or 



io + 3HAo - ^-^Ai = 



(140) 



Substituting (|140p in (jl38p . we have the dynamical field equation of Aj in 
an expanding universe as 



For the transverse direction QiAi = 0, we have 

Aj + hAj + + m^) Aj = 



(141) 



(142) 



where this is valid for photons (m = 0) and massive vector bosons in the 
transverse direction (A = 1,2). 
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For the parallel direction (A = 3,m 7^ only), we have 

-,2 , , J2 



It is impossible to solve equation p43p exactly [13]. However, at late time 
during inflation, a(t) grows more or less exponentially. Therefore, the second 
term proportional to Aj in the equation above may be negligible. Hence, 
the vector field can be written as 

A(x,t) = /"Aj;pq-e,(g,A)a(q,A)^g(t) 

+e-*^-e:(g,A)a*(q,A)^*(t)] (144) 

where 

2 

^ e*{q, X)ej{q, A') = 6ij - (145) 

A,A'=1 

for photons m = and 

3 

e:(g",A)e,(g,A')^% (146) 

A,A'=1 

for massive vector bosons during late time inflation. Therefore, Aq{t) is the 
solution that satisfies 

I (a(i)|A(i)) + W + m^Mt) = (147) 

To solve the equation above at a general momentum q, we can work in the 
conformal time r. Hence the massive gauge field equation in an inflating 
universe is 

2/1 , ^2 



' + U' + ^)A = (148) 



r 



where r = ^. 



We see from eq. ()148p that in the limit of m = 0, the solution of a 
massless vector field is a plane wave. This solution is the same as those of 
the conformal scalar and massless fermion. The positive mode solution for 
a massless vector field at general wavelength is 

Aq{t) = 1 e-''^\ m = (149) 

(27r)2^ 
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For a massive vector field m ^ 0, the field equation (|148p has the Bessel's 
equation type [11] 

// / 9 4i/^ — 1 , , , 

Therefore, the general solution of a massive vector field is 

Ag{T) = £g^Hi^\-qT)+J^gV^HP{-qT),m^ (151) 

where 

- - r2 (152) 
4 

Since we want the solution to match the positive solution at deep inside 
the horizon e"*"^"^, only H^\x) but not H^\x) gives an e"*'^'^ factor in the 
large limit. Hence, J^g = and 

Ag{t) = £gV^Hi'\-qT) (153) 

A normalized constant £g is chosen to match with the solution at deep inside 
the horizon. Inside the horizon, the positive frequency solution is the same 
as that in flat space, which is, 

A(*) ^ —j-^=exp(-i I ojg{T')dT'\ (154) 

J —oo 



where Wg(r) = \/ q"^ + (ma)'^. With the property of Hankel's function in the 
asymptotic limit, |x| ^ oo 

4^)(x)-/^exp(z(x-^-J)) (155) 



Since we now allow the existence of a mass term which can be either large 
or small, the normalized constants £g can be a function of mass and this 
may affect the result of the correlation function. 



During inflation, the mass of the matter could be large due to the in- 
teraction of matter with the inflaton (p. As mentioned earlier, the slow roll 
condition of some inflationary theories requires m = Mpi. This can 
make the mass term as large as Mpi and may affect the final result of the 
correlation function. To determine the mass dependent coefficient £g, we 
match the solution with that inside the horizon. From eqs. (I153p . (1154p . 
and (jl55p . we have the mass dependent coefficient £g as 



g(^) = ^^ex(i+2-) (156) 
2(27r)2 
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From eqs. (jl53p and (jl56p . we therefore have the massive mode solution of 
the gauge field that we will use for the propagator as 



Ait) = -^ex(i+2.)^^(i)(_^,) (157) 
2(27r)2 



where = y ^ — r^. 

To calculate what the interaction vertices are, we need to expand beyond 
quadratic order in fluctuations. As derived in the appendix, the cubic order 
is 

H^AAit) = -J d^'xeHa'^f^^A^ + ^{d:Aj - djAif)v-^C, (158) 

where we use the energy momentum tensor in [7] and choose gauge Aq = Q 
for m = 0. For m 7^ 0, we can solve the constraint equation of Aq and plug 
it back into the action. As seen from eq. ()137p . Aq is decaying as -^^^ after 
horizon exit. Therefore, we can approximate the massive vertices as that 
of the massless case. It is only the propagators that will be different from 
massless case. To continue, we use the general formula in eq. (|40p with the 
replacement 



- ^A^,+^,{{^^A,f-^,A,^,Ai) (159) 



Cg(ti,2) ^ -4(ii,2)/g' (160) 
V(t) = -eHa^{t) (161) 



Since the purely electric term (through (Ai)'^), purely magnetic term (through 
(diAj)^ - {diAj){djAi)), and two cross terms arise when we evaluate the 
commutator [Hi, [H2,Q]], eq. becomes 

■^Ap{h)Ap'{ti)A;{t2)Ap{t2) + 

-^Ap{h)Ap>ih)A;{t2)A;,{t2) + 



-^Ap{h)Ap>{h)A;{t2)A;,{t2) (162) 
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and C, remains the same as in eq. ()69p because we are calculating the same 
correlation function (CC) with various kinds of matter loops. 

We now need to calculate the polarization factor Vi for i = 1...4. 
The Vi factor comes from the purely electric field term which is 



A,A'=1 

= l + {p-p'f (163) 



The V2 factor comes from the purely magnetic field term which is 
2 



^2 = ^ |(P • p')(ep,A • ep/,A') - (P • ep',A')(p' • ep,A)|^ 

A,A'=1 

= p'^p\l + (p ■ pf) (164) 
The ^3^4 factors come from the cross terms which are 

V3 = V4 = -2p • p' (165) 
Substitute these into eq. (|30]) . we get 

J d=^xe*'i-(^-''')(c(x,t)C(x',t))^ =-S{2^f j (iVV<5^(q + P + p) 

X / dt2e2H2al j dtieiHialUeizMjiXim) 

J— 00 J —00 

where 

Ma = (l + {p-pfy^^a2^Ap{ti)Ap'{ti)A*p{t2)A*p>{t2) 

+pV'(i + {p-pf)a^%^Ap{h)Ap>{ti)A;{t2)A;,{t2) 
-2pp\p ■ p')a:[*a^'^Ap{ti)Ap'{ti)Al{t2)Al,{t2) 
-2pp'{p-p')a^^a'^^Ap{ti)Ap:{h)Al{t2)Al,{t2) (167) 

Since the solution of the massless vector field is just a plane wave, 

A,{t) = -^AS) (168) 
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eq. (|167p can be simplified further as 

= tK, , (i + (p-pO)V-(^+^')(-^-^) (169) 

Substituting the Z part in eq. ([77j) and gauge field part in eq. (jl69p into 
eq. (I166p . we have the correlation function due to massless vector fields: 

J d3xe^'i-(^-^')<|c(x,t)C(x,t))^ = -16{27rfe^\CX J 

xS^q + p + p')\A;\'\A;W{i + iP ■ p)yT (170) 

where 

r = Re f (ir2(e-^''^2 _gigr2)g*(p+p')r2 P ^^^g-%+P+p')n 

J — oo J — oo 

1 



(171) 



2q{q + p + p' 

and the constant coefficients after horizon exit are 

Substituting eqs. (|17ip and (|172p into (|170p . we have the loop power spec- 
trum due to a massless vector field as 

/ ^'--'-''•'(C(x.*)<(x'.*)),^^ = 1^ / W 



x53(q + p + pO /^^ \ 1 + {P-P')) (173) 

Notice that the first term is the same as that of a massless minimal coupled 
scalar loop. The p-p' terms come from summation over polarization vectors. 
We can follow the same method of dimensional regularization shown in the 
earlier section. We have the finite part of correlation function due to massless 
gauge field as 
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Notice that the result is smaller than the classical result by a factor of GH^ 
in order of magnitude. The numerical coefficient is slightly more than that 
of a scalar loop in [2] because there are additional polarization factors. 



VII. MASSIVE VECTOR BOSON 

We consider the late time mode solution of the massive vector field because 
the mode solution for the propagator is no longer a simple plane wave as in 
the massless case. For |x| = | — pt\ — > 0, we have 

By definition of Hankel's function, 

H'^^Xx) = -^(j-Ax) - e-^^'Ux)) (176) 
I sm VK \ J 

Hence, the late time behavior of mode solution approaches 

The exact solution in eq. ()157p approaches 

Aq{t) = Cga^+ + Vga^- (178) 

at late time where 



1 1 /l m2 ^ ^ 



and 



' 2(2^)IV^sini/^r(l-r/) V q J ' ' 

' 2(2^)i^/:H^sini/^r(i + z.) V g ; ^ > 

The time derivative of the solution contributes the same power of a as 

Aq{t) = H(x+Cga^+ + X^Vqa^-) (182) 
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We see from the equation above that the time derivative of the propagator 
in massive theories contributes an additional factor 



A, - A±/(t) ^ ( - 1 ± - (183) 

where f{t) has the same power of t as in Aq. Therefore, there will be a term 
like 

Ap{h)Ap>{h)Al{t2)ACp,{t2) ^ 0{Xi)F[tiM (184) 
in the loop. Since at large mass limit, 

Ai - ^ (185) 

the time derivative propagators give an additional factor of (A-t)^ for four 
fields. Because of eq. (|185p . the loop spectrum does not seem to get sup- 
pressed by an additional factor of G, but may get suppressed by the constant 
coefficient at large mass in eqs. (11801) and (I18ip or the results of the loop 
integrals. We therefore investigate whether there is true suppression or not. 

To analyze in more detail, we start from the interaction in eq. ()158p . 
Therefore, the massive vector field loop contributes 

MA,m^o = n*(t2)n(ti) (186) 



where 



n(i) = -^^Ai^p^x[t)Ai^p,^x'[t) + -^^{^p^^ 

-Pip[Aj^p,x{t)Aj^p,^x'{t)) (187) 



The exact solution of a massive vector field involves Hankel's functions which 
are rather complicated to integrate over time and momentums. However, we 
can get some ideas about what the momentum dependence of the observable 
spectrum is by considering the long wavelength mode solutions. We see from 
eqs. ()178p and (|182p that Aq gives the same power of a as Aq at late time. 
Therefore, we can keep the most leading order term as the universe rapidly 
expands 

Ma -4-^A*{t2)Ap*{t2)Ap{ti)Ap.{ti) (188) 

0-2^1 
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This means that the (massive) electric-hke term is more dominating than 
the magnetic- hke term after horizon exit. This result is different from the 
result of massless vector fields in which all electric and magnetic terms are 
equally important. Substituting eq. (jl88p into eq. (fl0|) . we have 

j d=^xe*'i-("^"')(c(x,t)C(x,t))^^^^^ -24(2^)9 y dVV<5'(q + P + P ) 

dt2e2H2al j dtieiHia\Re(zAp\t2)A/ {t2)Ap{ti)Ap'{ti)\im) 

where Z is the contribution from the part which is still the same as that 
in eq. ([77|) . Hence, 

J (iW'i-(^-^')^C(x,t)C(x',t))^ ^ -24(2^)9|e|2|Q°|4y" (fpd^p'5^{q + p + p')Xl90) 
where 



T 



Re 



dt2a2{e~ 



■igT2 _ „iqr2 



e'^^^)Ap {t2)Ap> {t2) 



t2 



dtiaie-''i^'Ap{ti)Ap>{ti] 



(191) 



With the time derivative of the late time mode solution in eq. (jl82p . we 
have 



Ap{ti)Ap'{ti) 



}?^_CpCpi a-^ ^ + X^DpTDpia^ 



+X+X-{CpT>p' + T>pCpi)a^ 



\++\- 



(192) 



Therefore, the ti integral is 

'■t2 



dtiaie-'''^^Ap{ti)Api{ti) = / dhe'''^^' XlCpCp^aj 



t2 



-.21' 



+XtVpVp>a^^'' + X+X^{CpVp> + VpCp> 



H 



2u 



2v 

c+a2 



Xl 
2^ 



c_a2 — X+X-CoEi{—iqT2) (193) 



for 2i/ 7^ and 



Co 
c+ 
c_ 



CplDpt ~\- TDpCp' 

^p^p' 
VpVpi 



(194) 
(195) 
(196) 
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and A± = — |zb(z/ = Y| — ^) which can be either real or complex or zero, 
depending on its mass when compared to the expansion rate H. 

Small Mass: m < ^ 



For m < ^ is real. Therefore, 

-^<A+<0,-l<A_<-^ 
and \Cq\^ and \'Dq^ are q-dependent, depending on its mass. Hence, 



(197) 



r 



H^Re 



A 



-c+a 



2u 



dt2{-2i) sin qT2 \\c\af + Alcloa + A+A-Cq 
A2 



-c_a 



-2v 



X+X-CoEi{-iqT2] 



(198) 



.21/ ' ^ 2u 

We can see that there is no contribution from the terms proportional to 
|c+p, and |c_p because they are all real. With the factor i in the integrand, 
the contribution is purely imaginary. Hence there is no contribution after 
taking the real part. Therefore, the terms that give non-zero result are 



T 



-Re / ismqT2 

J-oo ^2 



2u 



+ . -A+A-|cop)^z(-igr2) 



{-HT2f^ {-Ht2)- 



-2X+X^ilm{c*_c+) - + 



Xlclc. 



(199) 



{-HT2y- i-HT2)-^-. 

The integral above is still too complicated. However, when 2u < 1, 
singr2 — qT2 = — and Ei{—ix) ~ Inx when x — 0, the result of 

the time integral is less than or (Ina(t))^ at late time. Therefore, the 
dominant contribution comes from the first term of the integrand. Hence, 

(1 — rjlnqr) 



T ^ 2qHX-XXlm{cXcQ)- 



(200) 



for <C 2i/ < 1. Therefore, the correlation function due to a massive vector 
field loop is 

12(27r)-^A_A-^(87rGi/2)2^ 



j d3;^e^'i-(--')(c(x,i)C(x',t) 
(1 — In gr) 



q^an 

d^pd-^p'S^{ci + p + p')Im{clco) (201) 
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To calculate what Im{c\co) is, we see from eqs. ()180p and pSip that C*Pp 
is p-independent. Therefore, 



Im{c\co) 



where we use r(l — x)r(x) 



(2ff)2^r^(t/) [1,1 



p2u ~^ p/2i' 



and 



Im{C;Vp) 



1 



(202) 

(203) 
(204) 



Hence, 



3(2)2'^r2(i/)A_ A3 (87rGF2)2 (l + {2u - 1) In qr) 



loop,m<l^ 



2(2^)4i7l-2'^g5al-2'^ 

(f p(fp' 6^ {q + p + p') 



v(\-2vf 
1 



(205) 



To determine the momentum dependence q of the spectrum, we integrate 
over internal momentum circulated inside the loop. Following the simi- 
lar way as in the massive fermion section, we have the momentum dependent 
spectrum as 



24r^{u)\^Xl{8TTGH'^f (1-r/lngr) 

{2TT)^{2Ha{t))'i (l-r?)(2 + r?)7/2g3-r? 



where 



2u 



Am? 

1 TT — 1 



2m^ 



or 



2m2 1 



(206) 



(207) 



(208) 



We see that the departure from scale invariance is still small. 
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Critical Mass: m = ^ 



For m = u is zero. Therefore, 



A+ = A_ = -^ (209) 
and Cq and Dg are q- independent. Eq. (|193p becomes 

J ' dtiaie^*'?"Mp(ti)^'(ti) = (Cp + Pp) (Cp/ + Pp.)^i(-igr2)(210) 

Substituting equation above into eq. (|19ip . we therefore have 

r = -^|Cp + pp|2|Cp/ + pp/|2 (211) 

We see from eqs. (|180p and (jl8ip that the coefficients are all momentum 
independent when 1^ = 0. Since 

we need to use 1' Hospital's rule. From eqs. (jl80p and (jlSip . we have 



^ \ iyfrre 4 
lim Co + = o lim 

--0^ ^ 2(2^)t/ff--o 



sin UTTr{l — u) sin z^7rr(l + v) 



(213) 



Note that 



e^(^r e^(f(f)^ + .(f)'^-^) 



lim — = lim ^ , , ,^ , 

>o sm wKLyl — z/j i^— ►o tt cos z^7rr(l — i/j — sin V7ryj{l — z^)r(l — z^j 

i 



Similarly, 



2 



lim ' ^ = -- (215) 

i/^o sini/7rr(l + z^) 2 ^ ^ 



Therefore, 



lim (Cg + Pj = ^ , ^ (216) 
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We see that the coefficients are g- independent. Hence, eq. ()21ip becomes 



r 



vr 



1024(27r)6 

Therefore, the correlation function in eq. (jl90p becomes 



j d3a;e''i-(--')^C(x,t)C(x',t) 
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Large Mass: m > -f- 



For m > \± are complex conjugates of each other as 



(217) 



(218) 



Ai = A 



+ is, X- = A* 



IS 



(219) 



where u = is and s = y\jj^ — \\ is real. Notice that |Cgp and are 
q -independent. Therefore, eq. (jl93p becomes 



t2 



dtiaie"*^^Mp(ti)X'(*i) 



H 

2il 



-2is\X\'^ coEi{—iqT2) 



(220) 



We see that the time components of the first two terms are complex con- 
jugates of each other with different constant coefficients. Integrating over 
time t2 gives 



r = -—Re 



ij. ■ \*2 * —2is , \2 * 2is , |\|2 * 

dt2smqT2 X c+Og +A c_a2 +\X\ Cq 

2 



X^c+af' - X*^c-a2^'' - 2is\X\''coEi{-iqT2 



iu2 fO 



\X\^K 



dT2 . ,2 I |2 

smgr2 |c+| — |c_| 

00-^2 V 



2s\co\'^SiqT2^ (221) 



where the other terms vanish in the limit of t — > oo because of the oscillating 
behavior of the integrand e^*^*^. We can integrate further with the use of 
Mathematica 



J —c 



dT2 . ^. IT'' 

smqT2SiqT2 = — — 

oo 7-2 8 

° dT2 vr 
smgr2 = — 

oo T2 



J-oo T2 



(222) 
(223) 
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Therefore, eq. (|22ip becomes 

7r|A|4if2 



r 



4s 



2(|c+ 



2 I |2\ I I |2 

p - |c_p) + 7rs\co\ 



We now need to calculate what |c-|-^oP From eq. (|194|) . we have 



|c+| 



|col 



CI 2 1 1 2 
I p\ \ ^p' I 

|Pp|2|Pp,|2 



From eqs. (I180p and ()18ip . we have 



and 



TT 



4(27r)3i7|sin7rz^|2 |r(l - 



TT e 
4(27r)3i7|sini/7r|2 |r(l + i/)|2 



(224) 



(225) 
(226) 
(227) 



(228) 



(229) 



We can see that \Cg\'^ and \'Dq\'^ are momentum independent because u = is 
is purely imaginary in the large mass limit when m > Therefore, 



I |2 I |2 

|c+| — |c_| 



-2ns 



) 



16(27r)6iJ2| sinz^7r|4|r(l + 



(230) 



Using |r(l + is) 



|r(i 



sinh TVS 



and sin is = zsinhs for real s, the 



equation above is simplified as 



I |2 I |2 
C+ — C_ 



coth vrs 



4(27r)6s2/72 
where we use sinh2x = 2sinh j;coshx. Also, 



(231) 



|co| 



cos2(sln^; 



4(27r)6s2i?2 ginh2 vrs 



(232) 



Notice that the coefficients |c+|2 — |c_|2 are completely momentum indepen- 
dent and I Cop is nearly momentum independent (In j^^^ — > when (7 — > 0). 
Substituting eqs. (I23T]) . and (f232ll into eq. (f224ll . we have 



r 



8s(27r)6 L 



coth vrs vrcos2(sln£:; 
^ — + 



2s sinh vrs 



(233) 
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From eqs. ()190p and (|233p . we have 



J d3^e^'i-(^-^'){|c(x,t)C(x',t) 



loop 



37r(87rGg^)^|A| 
32s(27r)3g6 



X(53(q + p + p') 



p 



cothvrs 7rcos^(sln|:) 
1- 



2ssinh tts 



(234) 



The result of the momentum integrals p, p' gives the momentum dependence 
as g^, which cancels with the factor in eq. ()234p . We therefore have the 
approximated scale invariant spectrum after horizon exit as 



J d3xe^'i-(--')(c(x,t)C(x',t) 



8s3(27r)3g3 



loap,rn>- 



■ COth TTS (235) 



where X = + is and s 



1 1 



VIII. CONFORMAL SCALAR LOOP, INFLATON, AND 

GRAVITY 

We have learned that the spectrums of massless minimal coupled scalar, 
massless fermion, and massless vector fields loops all go as {8TTGH'^)'^q~^ In q. 
We would like to investigate whether this is also true for conformal scalar 
loop. The full action considered during inflation is 



1 



2 



1 



>C = --V^ + g'^'d^^d.ip + 2V{ip) + g'^'-d^xduX - iRx' (236) 



87rG 



where ip is an inflaton , x is additional conformal scalar matter in which 
(x) = 0) ^-iid C = g)0 for conformal and minimal couplings respectively. 
We consider this to see how the conformal scalar affects the spectrum {QC,) 
through the interaction with gravitational fluctuation. To arrive at the field 
equation of the conformal scalar field, we need the action up to the second 
order in the field fluctuations which is 

4'^ = y(x^-^-124^V) (237) 

where R = —12H^. Varying the second order of the action with respect to 
X, we have the field equation of the conformal scalar field as 

2 

Xg + 3Hx, + + 12^^/') Xg = (238) 
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Notice that if there is no extra term 12^H^x {^^ C = 0)) this is just a 
minimal coupled massless scalar, in which the dominant solution approaches 
a constant at late time. It is known that the massless minimal coupled scalar 
produces a scale free spectrum. We would like to investigate the momentum 
dependence of the power spectrum due to conformal scalar loops here when 
^ = |. Eq. (j238p can be solved exactly by re-scaling the field x = u/a. 
Hence, for the conformal scalar. 



/ a" 

Uq + - — + 2i?^a^ )ug = (239) 



During inflation, a ~ "7^7; therefore, the last two terms of eq. (j239p are 
cancelled. We arrive at a simple field equation of the conformal scalar field 

< + q^Uq = (240) 

Therefore, the solution to the equation above is just a simple plane wave 
valid to all wavelengths 

Xq{t) = uq{t)/a{t) = \ e"-'?- (241) 

where we choose the constant coefficients to match with the positive mode 
solution inside the horizon. From eq. (j24ip . the conformal scalar field 
correlation function to leading order is 

X(x,t)x(x',t)) = I d^qe^^<^-^'^\xM' 

„iq-(x-x') 1 



(27r)^ 2qa^{t) ^ ' 

We see that its momentum dependence is far from scale invariant at the 
classical level. However, we never observe the product of the scalar fields in 
CMB anisotropy but rather the correlation function of the temperature or 
density fluctuations which is related to the conserved quantity C,. Therefore, 
we study how the conformal scalar field affects the observable power spec- 
trum {QC) via the gravitational interactions at the quantum level. We can 
calculate the trilinear vertices due to the conformal scalar x a-nd gravity C,. 
They are 

^Cxx = -lC{drX? + adi[jl-^Ha^'^-^C)xdiX 

-a'{m'C + HC-^-§)x' (243) 
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where the last hne above contains the additional terms from the massless 
minimal coupled scalar. Those terms in the last line arise from the confor- 
mal term which is ^^^^-Ry^ . We see that the interactions above are rather 
complicated. As derived in the appendix, the interaction vertices can be 
written in a more compact form as 

Hc^^^it) = - J d?xeHa^{T^^ + a^T")V-^C (244) 

where T^^^ is the energy momentum tensor at the second order of arbitrary 
matter. The combination of time and space components of the energy mo- 
mentum tensor is [14] 



^UU ^^2J,^^ ^ 2(1 - S^x' + 12^ ^'x' + ^X;i " X;00 + 

-ex' (i?oo + ^-R + ^iR) (245) 

We can check that T^^ + a?T^^ = 2%^ for minimal coupling ^ = 0. For 
conformal coupling ^ = |, we have 

TOO + a'T^ = x' + \ - 2xx - H\') (246) 

where the conformal field equation (j238p is used and 

Roo + ^-R + 3CR = (247) 

due to ^00 = ^H^,Rii = -9a^H^ and R = -12H^ in de-Sitter phase 
inflation. Therefore, the trilinear interaction Hamiltonian of the conformal 
scalar x and gravity C is 

Hc^xit) = - j d'xeHa^ [x' + i(^^ - 2xx - H\')] V'H (248) 

To calculate the loop spectrum in the commutator [Hi, [H2, Q]], we can use 
the general formula in eq. (|40p with the replacement 

- x' + ^(^-2xX-i^V) (249) 

Cq{h,2) ^ -Uti,2)/q^ (250) 
V{t) = -eHa^{t) (251) 
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Hence, 



Tr*{t2)TT{ti) 



where 



■ ■ 1 •■ 1 •• 1 fPP' , rr2\ 

= Xp'Xp - ^Xp'Xp - ^Xp'Xp - ^[-^ + ^ jXp'Xp 



With 

Xgit) = 

many terms are cancelled. Therefore 



a / 



a 



Xq 



= ^[P^+ P'^ - "^PP'l Xpit)Xp'{t) 



(252) 

(253) 

(254) 
(255) 

(256) 



Hence, 



36(27r)6afafpp' 
From eq. (j40p . we have 



1 / \ 2 

^[P^ +P'^ - ^PP) Xl{t2)X*p'{t2)Xp{tl)Xp'{tl) 

\ , , + _ 4^y\2^»»(p+p')(n-r2) ^257) 



9af 02 



y (i3xe*'i-(^-''')^C(x,t)C(x',t))^ =-8{27rf J (iVV<5^(q + P + P ) 

dt2€2H2al I ' (itieii?iafi2e('z7r^7riV258) 

D J — CO 



where remains the same as in eq. ()77p . Substituting the 2^ part eq. (j77p 
and matter part eq. ()256p into eq. ()258p . we have the correlation function 
due to the conformal scalar field loop 



loop 
2 



2{2Tr)^{8TTGH' 
9^ 



2\2 



x5^(q + p + p') + p" - W Re / dt2a2{e-''^^' - e'^-^)xl{t2)xl'{t2) 



t2 



X / dhaie-'^^'xpiti)Xpiti) 



(259) 
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With the mode solution in eq. (j24ip and the result of the time integrations, 
we have the loop correlation function as 



X / d'^pdV'^ (q + p + p 



pp'{p + j/ + q) 



(260) 



With the dimensional regularization as done before, we have the conformal 
scalar loop correlation function as 



d3xe^'i-(''-^')/C(x,f)C(x',t) 



loop 



90(27r)3g3 



Inq + C 



(261) 



We see that it is nearly scale invariance and smaller than the classical result 
by a factor of ^"^^^/i^^"^ ■ 



IX. SUMMARY OF ALL RESULTS 

We study quantum effects of cosmological correlations due to the interactions 
of gravitational and matter fluctuations. It is shown that departures from 
scale invariance are never large, regardless of what kind of theories, what 
kind of matter, or what kind of inflaton potential V{ip) are used. 
The results in this paper may be compared with the Weinberg's result [2], 
Minimal Coupled Scalar Field Loops 



(CO 



m=0 



7r(87rG)^ff(tg) ^ - 
15(27r)3g3 



Ino + C 



(262) 



where H{tq) oc g is the expansion rate at the time of horizon exit. 



The results in this paper are 



Dirac Field Loops 



(CC)my=0 



4^(87rG)2//(tg)4 



15(27r)3g3 



hiq + C 



(263) 



24g3cosh2^ 



(264) 
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Gauge Field Loops 

(CC)m„=0 



4 r 



Inq + C 



(CO 



30(27r)3g3 

2AT'^{v)\_\\{STTGH{tqff (1 - 7?lngr) 

(27r)3[2i7(tg)a(t)]'' (1 - 7?)(2 + 7?)r/2g3-r, 



(265) 
(266) 



where A± = -i±z^, 0«;2z^ = Y^l-^<landr/ = l-2i/< 0.5. 

7r2(87rG)2/7(f,)4 



(CO, 



(CO^„>f 



where A = — | + is and s 



" 2 
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3^2(87rG)2/i'(tg)4|A|4 



^s3(27r) 



3„3 



COth TTS 



(267) 
(268) 



1 1 



IP 



2,Tl'^{MGfMpiH{tqf 



8(27r) 



3^3 



Conformal Scalar Field Loops 

■K{ST:GfH{tqY 



(CO 



90(27r)3g3 



ln« + C 



(269) 



(270) 



We see that even when the mass is as large as Mpi, the one-loop result is 
still naturally smaller than the classical one. Therefore, no fine tuning is 
needed. The result above is still valid in the realistic and the general poten- 
tial V {^p^'iI)iIj, Afj_A^). The reason is that we choose the gauge in which the 
inflaton does not fluctuate [7] {Sip = 0). 

Even when the additional interactions of inflaton and matters arise, the 
results in eqs. (j264p and (|268p do not change but rather the masses are 
shifted by 



ruf 
ml 



rrif + - |v.=o 



(271) 
(272) 

(273) 
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The is because the mass shift ( ^J^lv=0) c)a^(Mm U^=o), which is a function 
of the unperturbed inflaton only, does not change much during inflation. We 
therefore can approximate the unperturbed inflaton at the time of horizon 
exit (p{t) ~ <f{tq)- Hence there is no additional consequence to the momen- 
tum dependence of loop spectrums. 

Therefore, the spectrums are nearly scale invariant even if we add inter- 
actions of arbitrary matter and the inflaton to the interactions of matter and 
gravity. These results imply that we and the things around us did not come 
from nothing or an unknown scalar field as in conventional beliefs. Rather 
it points to the fact that we originated from quantum fluctuations due to 
the interactions between gravity and various matters during the time of Big 
Bang inflation. 
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APPENDIX: HIGHER ORDER FLUCTUATIONS 

This appendix is to clarify and derive interactions of matter and gravita- 
tional fluctuations used in the loop calculations. The method of expansion 
and quantization shown by Weinberg [2] has a more compact form than 
the direct expansion of matter and gravitational fluctuations. Following his 
method, we can extend the calculation to other matter such as fermion, 
gauge, and conformal scalar fields without many diflficulties. We would like 
to show the calculation in detail for the general reader. 

In cosmological fluctuations, we generally expand the gravity and an 
inflaton around a time dependent background such that 

5^^(x,t) = g^^{t) + 6g^^{x,t) (274) 
ip{x,t) ='f{t) + 6^{x,t) (275) 
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When we add any other kinds of matter which have unbroken symmetries, 
they can be expanded as 



M{x,t) = + 6M{x,t) (276) 

where M represents any additional matter such as fermions, vector bosons, 
and confer mal scalar fields. The perturbation to the metric around an FRW 
background can always be placed in the form of 

ds'^ = -(1 + E)dt'^ + 2a{t)Fidtdx' + a^{t){{l + A)6ij + B^ij)dx'dx^ (277) 

where we only consider the scalar mode which is the subject of interest here. 
The gauge invariant observable quantity is defined as 

«'-T-^ (278) 

to linear order. We see that C is a gauge invariant quantity that relates to 
both matter and gravitational fluctuations. There is a need for us to learn 
how to quantize such theories with minimum complication. 

Since the inflaton and gravity are related through Einstein's equation, 
we have some choices in choosing a gauge. It is found to be more convenient 
to choose a gauge such that the inflaton does not fluctuate [7] . Therefore, we 
can write down all the components of the gravitational fluctuations Sg^^i^ in 
terms of a single variable ^ by solving Einstein's equation in the Maldacena 
gauge Sip = B = 0. Prom the gravitational field equations and the energy 
conservation equations [9], 

= A- HE (279) 
= HE + 2{3H'^ + H)E + a-'^V'^A-A-6HA + 2a-^HV'^F{280) 

1/7 Q 

= (EH) -3HHE -a-^HV^F + -HA (281) 

2dV ' 2 ^ ' 

Solving the equations above, we therefore have 

A = 2C,E='^,F = -^ + eaV-^C, (282) 
H aH 

where e = — By eliminating E and F yields a differential equation for 
A: 
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This is sometimes known as Mukhanov equation[l], in which A = 2(. We 
can write the metric and its fluctuations in terms of as 



500 = - (l + §) = ^i^' - (284) 

90i = + ea^ V-^C) = Ni (285) 

gij = aHij(l + 2C^ =hij (286) 
The determinant of the metric is 

^^ = NVh = a^{l + ^)e^'^ (287) 
The gravitational, infiaton, and matter actions in Sip = gauge are 



C = 



0^ + 2V{ip) + ^ +CM{x,^i^,A^,A^') (288) 



where >Cm(x? V'V') A^A^) are the additional types of matter such as the con- 
formal scalar, fermion, and vector bosons that do not have the background. 
The first three terms give vertices of purely (. We are presently interested 
in the interactions of matter and gravitational fluctuations in the last term 
{jCm term) because, in general, the matter loops are larger than the ( loops 
by a factor of SttG. Therefore, the time dependent tri-linear vertices of 
general matter are 

where T^'^ is the energy momentum tensor of arbitrary matter evaluated at 
quadratic order in fluctuations. With the Bianchi Identity, 

we have 

\^{a^T^^) + aaT'' + diT'^ = (291) 
a-^ at 

where we use r*Q = 3H, F^- = aaSij , Fq^ = f ?q = F*^- = for the unperturbed 
FRW metric. Integrating by parts in space and using the Bianchi Identity 
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eq. (|29T]1 . eq. ([289]) becomes 



H(MM{t) 



I' 



c 



.H 

2\-7~2 /- / ^Irpii' 



Ha-^ a at 



where the term a (T^^ is cancelled. With the Mukhanov equation 



C + [3i/ + -]c-— C = o 

e 



Eq. (|292p is simplified as 



H^MM{t) = Z{t)+Y{t) 



where 



Zit) 
Y{t) 



6rr00 



c 



(292) 

(293) 

(294) 

(295) 
(296) 



The term Y{t) can be removed by the field redefinition of C = exp(— fy)C exp(iy) 
as mentioned in [2] . To see more clearly, for any interaction Hamiltonian of 
the the form ()294p . Eq. ([8]) can be put in the form 



(Qit)) 



dt]\f 



dt 



JV_1... 



dti 



Af=0 

x{[Hi{h),[Hi{t2),...[Hi{tN),Q\t)]. 



t2 



(297) 



where 



Hi{t) 



Z(t) + y(t) + ie-n*)(^^e^m 



-iY{t) 



Z{t) + i[Y{t),Z{t)] + -[Yit),Yit)] + ... 



(298) 



and 



Q'{t)e 



-iY{t) 



Q'it) + z[Y{t),Q\t)] - i[y(t), [Y{t),Q'{m + 



(299) 
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As mentioned in [2] , the redefinition of the operators is necessary. It is only 
products of the redefined field operators whose expectation values may be 
expected to give results that converge at late times. The results contributed 
from Y{t) part only give a sum of powers of q with divergent coeflacients, 
but with no logarithmic singularity in q. Therefore, ( used in the sections 
IV- VIII is a new redefined variable, in which we can safely calculate the 
contribution from the Z(t) part only. The calculation in this way is more 
simplified than the direct expansion of the fluctuations. 
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